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We have performed a search for the decays B + — > J/ippA and B° — > J/4>PP in a data set of (88. 9± 
1.0) x 10 6 T(AS) decays collected by the BABAR experiment at the PEP-II e + e~ storage ring at 
the Stanford Linear Accelerator Center. Four charged B candidates have been observed with an 
expected background of 0.21 ± 0.14 events. The corresponding branching fraction is (12^g) X 10~ 6 , 
where statistical and systematic uncertainties have been combined. The result can be interpreted 
as a 90% confidence level (CL) upper limit of 26 x 10" 6 . We also find one B° candidate, with an 
expected background of 0.64 ± 0.17 events, implying a 90% CL upper limit of 1.9 x 10~ 6 . 

PACS numbers: 13.20.He, 12.39.Mk, 12.39.Jh 

Studies of the inclusive production of charmonium published by CLEO Q and BABAR 0, and preliminary 
mesons in B decays at the T{AS) resonance have been results have been presented by Belle |3| One of the in- 
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FIG. 1: Center-of-mass momentum of J/ip mesons produced 
directly in B decays (points). The histogram is the sum of the 
color-octet component from a recent NRQCD calculation 4] 
(dashed line), which includes multi-body final states, and the 
color-singlet J/ip component from simulation || (dotted 
line) . 



teresting features observed by all three collaborations is 
an excess of J/ip mesons at low momentum in the e + e~ 
center-of-mass frame, pcm, when compared to distribu- 
tions predicted by non-relativistic QCD calculations 0. 
Figure ^ (from Ref. 0) shows pcm for J/ip mesons pro- 
duced in B decay after subtraction of the component due 
to the decay of heavier charmonium states. The excess 
below 0.8GeV/c corresponds to a branching fraction of 
approximately 6 x 10 -4 , 8% of the total direct J/ip pro- 
duction. 

Possible sources of the excess include an intrinsic 
charm component of the B or the production of an 
sdg hybrid Q in conjunction with a J/tp . Another pos- 
sibility H is that the excess is from decays of the form 
B — * J/ip baryon anti-baryon. The rate of these decays 
could be enhanced by the intermediate production of an 
exotic state allowed by QCD but not yet observed, in- 
cluding nuclear-bound quarkonium (a cc pair bound to a 
nucleon), baryonium (a baryon- antibaryon bound state), 
or a pentaquark (a baryon containing five quarks). If 
such resonances were narrow, the other particle in the 
decay would be monoenergetic in the B rest frame. Note 
that the J/ip spectrum in Fig. ^ would not directly dis- 
play such narrow distributions because it is measured in 
the e + e~ center-of-mass frame. The difference between 
Pcm and p*, the J/ip momentum in the B rest frame, 
has an RMS of 0.12 GeV/c due to the motion of the B. 

This Letter presents searches for the decays 
B + — > J/ippA and B° — ► J/ippp in a sample of 81.9 fb _1 
collected by the BABAR detector. Note that the latter 
decay is Cabibbo suppressed relative to the former. 
Charge conjugation is implied throughout. 

BABAR operates at the PEP-II e + e~ storage ring, 
which collides 9.0 GeV electrons on 3.1 GeV positrons to 
create a center-of-mass system with energy 10.58 GeV 
moving along the z axis with a Lorentz boost of (3"f = 



0.55. T(45) production makes up approximately 23% of 
the total hadronic cross section. 

The BABAR detector is described in detail in Ref. [tj. 
The trajectories of charged particles are reconstructed 
and their momenta measured with two detector systems 
located in a 1.5-T solenoidal magnetic field: a five-layer, 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH). The tracking fiducial volume cov- 
ers the polar angular region 0.41 < 8 < 2.54 rad, which is 
86% of the solid angle in the center-of-mass frame. The 
transverse momentum resolution is 0.49% at 0.3GeV/c 
and 0.59% at 1 GeV/c. 

The energies deposited by charged tracks and photons 
are measured by a CsI(Tl) calorimeter (EMC) in the 
fiducial volume 0.41 < 6 < 2.41 rad (84% of the center- 
of-mass solid angle) with energy resolution at 1 GeV of 
2.6%. Muons arc detected in the IFR, a multilayer device 
of resistive plate chambers located in the flux return of 
the solenoid. The DIRC, a Cherenkov radiation detector, 
is used to identify charged particles. 

We select B candidates of interest in a BB-enriched 
sample. Events in the sample are required to have vis- 
ible energy E greater than 4.5 GeV and a ratio of the 
second to the zeroth Fox- Wolfram moment 0| > R2 > less 
than 0.5. Both E and R2 are calculated from tracks and 
neutral energy deposits in the respective fiducial volumes 
noted above. The same tracks are used to construct a pri- 
mary event vertex, which is required to be located within 
6 cm of the beam spot in z and within 0.5 cm of the beam 
line. The beam spot RMS size is approximately 0.9 cm 
in z, 120 /an horizontally, and 5.6 fim vertically. 

There must be at least three tracks in the fiducial vol- 
ume satisfying the following quality criteria: they must 
have transverse momentum greater than O.lGeV/c, mo- 
mentum less than 10 GeV/c, at least 12 hits in the DCH, 
and approach within 10 cm of the beam spot in z and 
within 1.5 cm of the beam line. 

Studies with simulated data indicate that these criteria 
are satisfied by 96% of generic BB events. 

B + — > J/ippA candidates are formed by combining 
J/ip, proton, and A candidates. J/ip candidates must 
have mass in the range 2.950-3.130 GeV/c 2 if recon- 
structed in the e + e~ final state or 3.060-3.130 GeV/c 2 
in /it + /i~ • 

One of the two electrons from the J/ip must satisfy 
the following ("tight") requirements. It must have an 
energy deposit in the EMC between 89% and 120% of 
its momentum, a Cherenkov angle in the DIRC within 
3er of expectation for an electron, a lateral moment of 
the energy deposit Hll. LAT. between 0.1 and 0.6, an 
A42 Zernike moment |l2j less than 0.11, and an energy 
loss in the DCH consistent with expectation. Less strin- 
gent ( "loose" ) requirements are imposed in the selection 
of the second electron: we require an energy deposit in 
the EMC of at least 65% of its momentum and place 
a less restrictive requirement on DCH energy, with no 
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requirements on LAT or A^- Whenever possible, pho- 
tons radiated by an electron traversing material prior to 
the DCH (0.04 radiation lengths at normal incidence) are 
combined with the track 

At 1.5GeV/c, a typical lepton momentum, the tighter 
criteria have an efficiency of 91% with a pion misidentifi- 
cation probability of 0.13%. The looser criteria give 98% 
efficiency with 3% pion misidentification. 

Muon candidates must deposit less than 0.5 GeV in the 
EMC (2.3 times the minimum-ionizing peak) and have a 
pattern of hits in the IFR consistent with the trajectory 
of a muon. The total amount of material penetrated 
must be greater than 2 interaction lengths and must be 
within 2 interaction lengths of the value expected for a 
muon. The muon identification efficiency at 1.5 GeV/c is 
77% with a pion misidentification probability of 11%. 

Proton candidates are selected with a likelihood 
method that uses the energy deposited in the SVT and 
the DCH, and the Cherenkov angle and number of pho- 
tons observed in the DIRC. They are also required to fail 
the tight electron identification criteria. At a typical mo- 
mentum of 300 MeV/c, the selection efficiency is greater 
than 98% with a kaon misidentification probability less 
than 1%. 

The A is reconstructed from a proton, which must sat- 
isfy the above criteria, and an oppositely charged track, 
assumed to be a pion. It must have mass between 1.10 
and 1.14GeV/c 2 , and a vertex that is separated from the 
J/ip vertex by at least 2 mm. The angle between the A 
momentum and the vector from the J/ip vertex to the A 
vertex must be less than 90° in the laboratory frame. 

Geometrical vertex fits are performed on the resulting 
B + candidates, of which approximately 68% are rejected 
by a requirement on the quality of the fit. 

£>° — > J/ippp candidates are formed from J/ip candi- 
dates and an oppositely-charged pair of proton candi- 
dates. Approximately 83% of resulting candidates fail a 
requirement on the quality of a vertex fit. 

We use two nearly-independent kinematic variables |9j 
to categorize B candidates: the difference between the 
reconstructed and expected energy of the B candidate 
in the e + e~ center-of-mass frame, AE ~ (qr ■ qB — 
s/2)/^/i, and the beam-energy substituted mass, rriES = 

(0.5s + pb ■ pr) 2 /Ej. — p 2 B . The four-momentum of 
the e + e~ initial state, obtained from the beam momenta, 
is qr = {EtiPt), and s = \qr\ 2 - The four-momentum of 
the reconstructed B candidate, qB = (Eb,Pb), is found 
by summing the four-momenta of the three daughters, 
with daughter masses constrained to accepted values . 

The "analysis window" AW is defined by 5.2 < toes < 
5.3GeV/c 2 and -0.10 < AE < 0.25 GeV (B+ candi- 
dates) and -0.25 < AE < 0.25 GeV (B° candidates). 
The AE range is smaller for the charged candidates due 
to a kinematic cutoff in the B + — > J/ippA decay. Only 
candidates in the AW are considered in the analysis. Ap- 
proximately 15% of B + events and 1.5% of B° events 



contain more than one candidate, in which case we select 
the one with the lowest \AE\. 

For signal events, (AE) « and (to-es) ~ Mb. We de- 
fine a signal ellipse by [(to E s — M B ) /a m ) 2 + [AE/cte} 2 < 
S 2 , where the resolutions a m and <je are estimated from 
simulated data to be 3.1MeV/c 2 and 6.5 MeV, respec- 
tively, for B+ -> J/ippA, and 2.7MeV/c 2 _ and 5.5 MeV for 
B a -> J/ippp. S = 2.4 for B+ -> J/0 pA and S = 2.2 for 
B° -> J/ippp. 

The selection criteria for charged and neutral B can- 
didates, including the values for S, have been chosen to 
minimize the 90% CL upper limit expected in the ab- 
sence of real signal, based on simulated signal and back- 
ground events. Approximately 90% of the background 
events satisfying the criteria are combinatorial BB, in 
which tracks from the decays of both B mesons are used 
to form the candidate. The rest are continuum (non-BB) 
events. Both components are distributed throughout the 
AW, and neither peaks in the the signal of either AE or 

TOES- _ 

We use simulated B + — > J/ippA and B° — » J/ippp 
events to measure the selection efficiency. The simula- 
tion does not include exotic QCD bound states. We study 
the accuracy of the simulation of the detector response 
by comparing data and simulated background events in 
samples similar to the final selection. We compare the 
number of J/0 mesons reconstructed in B° —> J/ippp can- 
didates in which only one proton satisfies the identifica- 
tion criteria, and we compare the number of A baryons 
reconstructed in B + — ► J/ippA candidates in which the 
proton daughter of the B + is required to fail the cri- 
teria. Based on these studies, we apply multiplicative 
corrections to the efficiency of 0.97 ± 0.06 for J/0 recon- 
struction and 0.86 ± 0.14 for A reconstruction. We also 
compare the distributions of the % 2 of the B vertex for 
candidates satisfying all other criteria and obtain correc- 
tions of 0.98 ± 0.02 for B+ -> J/iPpA and 0.90 ± 0.10 for 
B° -> J/ippp. _ 

The efficiency for B + — > J/ippA, with the J/0 decaying 
to e + e~ or [i + fi~ and A decaying to p7r + , is 0.049 ± 0.009. 
The 18% fractional uncertainty includes 16% from A re- 
construction, 6% from the J/0, 3% from statistical un- 
certainty in the simulation, 2% from the % 2 correction, 
and 1% uncertainty on proton reconstruction efficiency. 
Approximately 25% of signal events satisfying all other 
criteria are reconstructed outside the signal ellipse. 

The efficiency for B° — ► J/ippp with the J/ip decaying 
to e+e- or /j,+fj,- is 0.184 ± 0.024. The 13% uncertainty 
includes 6% from J/ip reconstruction, 2% for statistical 
uncertainty in the simulation, 11% for the % 2 correction, 
and 3% for proton reconstruction. 

We use world average values 01 f° r B(J/ip — > e + e~), 
B(J/ip -> and B(A -> pn~). 

We estimate the mean expected background in the sig- 
nal ellipse (hb) from the number Na elsewhere in the 
AW: [is = f • Na- We obtain /, the proportionality 
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FIG. 2: (a) Distribution of B + —> J/ippA candidates in the 
AE-m,ES plane, with the signal ellipse and its projection 
in each dimension (dashed lines). Histogram of candidates 
within marked bands in (b) jbes and (c) AE. Plots (d)-(f) 
show similar quantities for B° — » J/ip pp. 



constant, from a larger sample in which only one proton 
satisfies the proton identification criteria. We perform 
a Kolmogorov test [3] to verify that the distribution of 
candidates in the AE'-mEs plane is similar to the stan- 
dard selection. Comparing the regions outside the ellipse, 
the test gives a probability of 0.52 for B + — > J/tppA and 
0.36 for B° -> J/ippp. We obtain / = 0.0054 ± 0.0035 
(B + ) and / = 0.0051 ± 0.0013 (B°). The uncertainties 
are largely statistical, but include a component (16% for 
B + and 2% for B°) due to differences in the number of 
events with multiple candidates. 

For B + — > J/ippA, Na = 39, implying an expected 
background of 0.21 ± 0.14 events. We observe four can- 
didates in the signal ellipse (Fig. EJ. The probability 
of observing > 4 candidates when expecting 0.21 ± 0.14 
is 2.5 x 10~ 4 . Three of the four are positively charged. 
Two of the four J/ip mesons decay to e + e~ and two to 

To interpret this result as a B + branching fraction B, 
we undertake a Bayesian analysis with a uniform prior 
above zero. We define the likelihood for B as the prob- 
ability of observing exactly four events, including un- 
certainties on the expected background, signal efficiency, 
secondary branching fractions, and number of Y(AS) de- 
cays, (88.9 ± 1.0) x 10 6 . We assume the branching frac- 
tions B(Y(AS) -> B+B~) = B(Y(AS) -» B°B°) = 0.5. 

The central value for B is the peak of the likelihood 
function. We obtain "±lcr" uncertainties from a con- 
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FIG. 3: Momentum in the B rest frame of the (a) J/xj) , 
(b) proton, and (c) A daughters of the four B + — > J/ippA 
candidates. 



fidence interval that encloses 68.3% of the area of the 
likelihood function, selected such that the likelihoods 
for all values of B in the interval are larger than the 
likelihoods outside. The result is B(B + — * J/ippA) = 
(11.6±|-J) x 10" 6 . We similarly obtain a 90% CL up- 
per limit of 26 x 10~ 6 . 

If we consider only the statistical uncertainty, the re- 
sult would be B(B+ -> J/ippA) = (11.6l^) x 10~ 6 . Sub- 
tracting these uncertainties in quadrature would indicate 
contributions from systematic errors of 4.2 x 10~ 6 and 
1.8 x 10~ 6 on the upper and lower sides respectively. The 
systematic error arises almost entirely from the uncer- 
tainty on the signal efficiency. 

The creation of a narrow QCD exotic bound state as 
an intermediate resonance in the B + decay would be re- 
flected as a narrow p* distribution of the other decay 
daughter. We do not observe any significant clustering 
in the p* distributions of the J/tp , proton, or A daughters 
of the four B + candidates (Fig.|2J). The resolution in p* 
is a - 20 MeV/c. 

For B° — » J/tppp, there are 126 events outside the 
signal ellipse, indicating an expected background of 
0.64 ±0.17 events, and one event in the ellipse. Fol- 
lowing the procedure described for B + — > J/tppA, and 
again assuming a uniform prior above 0, we obtain 
B(B° -► J/tppp) < 1.9 x 10~ 6 (90% CL). This limit is 
dominated by statistical uncertainty. 

In summary, we observe four B + — > J/ippA candi- 
dates in a data set of (88.9 ± 1.0) x 10 6 T(45) decays. 
The probability of the expected charged B background, 
0.21 ± 0.14 events, producing > 4 events is 2.5 x 10~ 4 . 
The branching fraction is (12^g) x 10~ 6 , where the un- 
certainty includes both statistical and systematic compo- 
nents. This result can be interpreted as a 90% CL upper 
limit of 26 x 10" 6 . 

We observe one B° — > J/ipPP candidate with an ex- 
pected background of 0.64 ± 0.17, and determine a 90% 
CL upper limit of 1.9 x 10~ 6 on the branching fraction. 

Neither final state makes a significant contribution to 
the observed excess of J/ip mesons in inclusive B decay. 
The momentum distributions of the B + daughters do 
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not provide evidence for QCD exotic particles produced 
as narrow intermediate states. 
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